A chiral spin liquid wave function and the Lieb-Schulz-Mattis theorem 
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We study a chiral spin liquid wave function denned as a Gutwziller projected BCS state with 
a complex pairing function. After projection, spontaneous dimerization is found for any odd but 
finite number of chains, thus satisfying the Lieb-Schultz-Mattis theorem, whereas for even number 
of chains there is no dimerization. The two-dimensional thermodynamic limit is consistently reached 
for large number of chains since the dimer order parameter vanishes in this limit. This property 
clearly supports the possibility of a spin liquid ground state in two dimensions with a gap to all 
physical excitations and with no broken translation symmetry. 
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A long time after its first proposal, [1| the existence of 
a spin liquid ground state (GS) in two-dimensional (2D) 
quantum spin one-half models is still a very controversial 
issue. This is mainly because all one-dimensional (ID) or 
quasi-lD spin models that can be solved exactly, cither 
numerically or analytically, 0, El display a gap to the 
spin excitations only when a broken translation symme- 
try is found in GS (e.g., in spin-Peierls systems or 
when the unit cell contains an even number of spin 1/2 
electrons (e.g., in the two-chain Heisenberg model |j]). 
Hence, in these models the electronic correlations do not 
play a crucial role since their GS can be adiabatically con- 
nected to a band insulator without any transition. This 
important property of insulators, which clearly holds in 
ID systems, has been speculated to be generally valid 
even in higher dimensions, as it appears to follow from a 
general result, the Lieb-Schulz-Mattis (LSM) theorem, |f| 
whose range of validity has been extended to more inter- 
esting 2D cases. 0,0 

Recently, there has been an intense theoretical and nu- 
merical investigation of non-magnetic wave functions ob- 
tained after Gutzwiller projection of the GS of a BCS 
Hamiltonian. [a, l2l llfll lUj O n a rectangular L x x L y lat- 
tice, this state can be written in the following general 
form: 
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where N is the number of electrons (equal to the number 
of sites, i.e., N = L x x L y ). Pq is the Gutzwiller pro- 
jection onto the subspace of no doubly occupied sites, 
and c\. | and cj. . are creation operators of a spin up or a 
spin down electron, respectively. These are defined in a 
plane-wave state with momentum fc allowed by the cho- 
sen boundary conditions: periodic (PBC) or antiperiodic 
(APBC) in each direction. It is worth noting that, af- 
ter the projection, the wave function Q corresponds to 
PBC on the spin Hamiltonian, regardless the choice of 
boundary conditions on the electronic states. The pair- 



ing function fk can be easily related to the gap function 
Afe and the bare dispersion ej, = — 2(cosfc x + cosk y ) of 
the BCS Hamiltonian: 
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by means of the simple relation = Ak/(ek+Ek), where 
the gap function A*, can be in general any complex func- 
tion even under inversion (At, = A_^), as required here 
for a singlet wave function, [9j and — y/\Ak\ 2 + e'f, 
represents the spin-half excitation energies of Hbcs- 

As clearly pointed out by Wen, ^3] in presence of a 
finite gap Abcs i n the thermodynamic limit, such that 
Ek > Abcs > 0, the corresponding BCS finite corre- 
lation length is expected to be robust under Gutzwiller 
projection. Here, we restrict to this class of non-magnetic 
states, considering the projected BCS (p-BCS) state 
which is obtained by a d + id gap function of the fol- 
lowing form: 



A x 2_ y 2 (cos k x — cos k y ) 
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This wave function breaks the time reversal symmetry 
T and the parity symmetry (x <-> y in 2D) P, whereas 
T <S> P is instead a well defined symmetry. Hence, this 
spin liquid wave function may have a non- vanishing value 
of the so-called chiral order parameter: ^jj 
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with d x = (1,0) and d y = (0, 1). Chiral spin liquids 
were introduced a long time ago, |l3lll4l| however, to our 
knowledge, this is the first attempt to represent this class 
of states in the p-BCS framework. On finite size systems, 
we take the real part of the complex wave function , so 
that all the finite-size symmetries, including parity, are 
satisfied and a spontaneously broken lattice symmetry 
can occur only in the thermodynamic limit. 
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In the following we consider in more detail the rela- 
tion of the p-BCS wave functions with the LSM the- 
orem. Given a short-range spin Hamiltonian H, on a 
L x x L y rectangle (PBC on the x direction are assumed) 
and a variational state l^o) with given momentum, we 
can define another variational state, = Olsm\^o), 
by means of the LSM operator: 
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where r = (x, y) indicates the position of each site on the 
lattice. The new variational state |-0q) has the following 
properties: [(| 

1. Its energy expectation value differs at most by 
0(L y /L x ) from the variational energy of \ipo). 

2. If L y is odd, regardless of the boundary conditions 
on the y direction, the momenta parallel to x cor- 
responding to \ipo) an d IV'o) differ by tt. Hence: 
(Wo) = 0. 

For ID or quasi-lD system with odd number of chains 
L y and vanishing aspect ratio (L y /L x — > for L x — > oo), 
by applying the LSM operator to the actual GS, it is pos- 
sible to construct an excitation of the system with mo- 
mentum (tt, 0) which becomes degenerate with the GS in 
the thermodynamic limit. This implies in turn either a 
gapless spectrum or, in presence of a finite gap, a two- 
fold degenerate GS with a doubling of the unit cell and 
a spontaneously broken translation symmetry. For in- 
stance, the presence of a singlet zero-energy excitation 
with momentum (tt, 0) is just a characteristic of sponta- 
neous spin-Peierls dimerization, as it appears for example 
in the Majumdar-Gosh chain. [2( This result, holding rig- 
orously in the limit of vanishing aspect ratio, has been 
argued to apply in general for 2D systems. In the 
following, we will show, with an explicit example, that 
this result in 2D does not necessarily imply spontaneous 
dimerization, but topological degeneracy of the GS. 

It is simple to show that 6 LS m\p~BCS) = |p-BCS'), 
namely the same type of wave function of Eq. JQ| is ob- 
tained, with the changes below: 
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where the new quantized momenta k = k + (tt/L x , 0) are 
obtained by interchanging PBC with APBC in the x di- 
rection only and Eq. JJJ) means that the pairing function 
is calculated with the old momenta k: fkc\ ^c_^ ^. By 

definition, the wave function \p— BCS') has therefore the 
same quantum numbers predicted by the LSM theorem, 
the change of momentum being implied by Eq. 0. The 
reason why the momentum of the wavefunction can 
be non zero for odd number of chains is indeed rather 
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FIG. 1: Dimer-dimer correlation functions as a function of 
the distance for three and four chains with PBC on both di- 
rections, using the wave function Q with A x 2_ y 2 = 2 and 



subtle but easy to verify. Indeed the x— translation op- 
erator with APBC translates all creation operators, but 
the ones belonging to the boundary are also multiplied 
by a phase factor (—1). This translation operator always 
leaves invariant the \p— BCS') wave function. However, 
for a spin state with one electron per site each configura- 
tion has always L y electrons at the boundary, so that the 
physical spin translation operation (defined with PBC), 
differs from the APBC one for an overall phase (— l) Ly , 
namely a momentum (tt, 0) for odd number of chains. 
Analogously, the excitations obtained by modifying only 
the boundary conditions in the BCS Hamiltonian (in the 
x and/or y direction), namely using Eq. © (and/or its 
equivalent for the y direction) and fk — ► /g, may display 
in 2D the topological degeneracy of this spin liquid wave 
function. OH [U 

We now assume that the the p-BCS wave function 
with the Afe given by Eq. J3Jl represents the GS of some 
short-range Hamiltonian. However, here we do not ad- 
dress the question whether this chiral wave function can 
be stabilized in some physical Hamiltonian. Certainly ex- 
plicit Hamiltonians with short-range off-diagonal matrix 
elements can be constructed by a simple inversion prob- 
lem scheme. Having a finite gap A^cs, this wave 
function describes a spin system with a finite correlation 
length, and consequently a finite triplet gap. The LSM 
theorem can be then applied in the geometries where it 
holds, like, for instance, the three-leg ladder with PBC 
in both directions. As shown in Fig. 0}, it is clear that 
spontaneous dimerization is obtained in the thermody- 
namic system for this geometry, as the dimer-dimer cor- 
relation functions on each chain, 

A(r — r ) = (S* Sr+dx^r' Sr'+dJ ; (8) 

behaves for large distance as (— \)( x ~ x )(9g p/3 6 + const., 
being Osp the dimer order parameter. [3|, [lfj In strong 
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FIG. 2: Left: Finite-size scaling of the dimer order parameter 
in the rectangular geometries with an odd number of chains. 
Right: Semi- log plot of the dimer order parameter as a func- 
tion of the number of chains. In both plots, the data refer 
to the chiral wave function with A x 2_ y 2 = 2 and A xy — 1, 
whereas full dots in the right panel refer to the short-range 
RVB (the latter wave function having £ — > oo in 2D). 
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FIG. 3: Chiral order parameter as a function of the number 
of chains. Its value is determined by studying the chiral- 
chiral correlation functions on the smallest parallel triangles 
at the maximum distance along the long direction. L x = 50 
is indeed large enough for converged results at fixed number 
of chains, as the wave function has a large gap Abcs — 1- 
The variational parameters are A x 2_ y 2 = 2 and A xy = 1. 



analogy with the one-dimensional Hciscnbcrg chain in 
the gapped phase, the broken translation symmetry al- 
lows the system to satisfy the LSM theorem. In fact, 
it implies two degenerate singlet states with momentum 
differing by (tt, 0). In contrast, on any even- leg ladder, 
where the LSM theorem does not imply the degeneracy, 
the p-BCS state does not break translational invariance, 
as illustrated in Fig. ^ for the four-leg system. Despite 
the dichotomy between the odd and even chain cases, the 
2D thermodynamic limit can be still consistently defined. 
In fact, as it is clearly shown in Fig. [21 though a finite 
dimer order parameter is obtained for any odd chain lad- 
der, the order parameter is exponentially decreasing with 
the number of chains [see Fig.|2j. This implies that the 
broken symmetry, which is correctly obtained for odd but 
finite number of chains, represents an irrelevant effect in 
the 2D thermodynamic system. Nonetheless, in the 2D 
system, the GS can possess degenerate topological exci- 
tations. For instance, the matrix element of the dimer 
operator with momentum (tt, 0) between the two dege- 
narate singlet states, which is finite on any finite number 
of chains, decreases exponentially with increasing L y , as 
it is bounded by order parameter (see Fig.[5J|. Remark- 
ably, the chiral order parameter remains instead a gen- 
uine feature of this variational wave function even in the 
2D limit, as shown in Fig. |21 where an order parameter 

Oc = \J (^c) — 0-03 was found for this wave function. 

We have given here a clear example that a spin liquid 
GS can be stable in 2D, and yet satisfying all the known 
constraints given by the LSM theorem. Indeed, sponta- 
neous broken translation symmetry is obtained for any 
odd number of chains, a remarkable feature since before 
projection the wave function is translational invariant. 
As the number of chains increases the chiral spin liquid 



appears only in 2D systems, where the spin-Pcicrls dimer 
order parameter converges to zero, and no broken trans- 
lation symmetry is implied in the thermodynamic limit. 

The chiral spin-liquid described by the p-BCS wave 
function is also consistent with a recent extension of 
the LSM to 2D systems with finite aspect ratio L y /L x . 
In this pointed out by Oshikawa, [7j the state 

|V»o) = Olsm\iPo}, is not necessarily degenerate with 
the starting wave function IV'o)) as in the usual LSM 
construction. However, through a well-defined adiabatic 
evolution - in analogy with Laughlin's treatment of the 
quantum Hall effect - one can obtain a different state 
|Vo)- with the same spatial quantum numbers of |Vo) 
and degenerate with IV'o)- If IV'o) is described by a p- 
BCS state {Q, it is easy, by a small change of the pairing 
function of the state \ip' ), i.e., — ► f%, to define a state 
IV'o) with the same momentum implied by LSM theorem 
but expected to be degenerate with IV'o)- In fact, the 
wave function \i/j' } can be obtained with the same BCS 
Hamiltonian with APBC in the x direction. Then, 
in presence of a finite gap in the excitation spectrum 
Abcs > 0, the wave functions |Vo) and |Vo) have the 
same value on any physical operator, and in particular 
the Hamiltonian, so that the considered states are de- 
generate in the thermodynamic limit. This is shown, for 
instance in Fig.^Ja) for the nearest-neighbor total energy 
contribution. 

The state |Vo)j obtained by the adiabatic evolution of 
the LSM excitation, \tp' ) = OlsmIVo), is no longer con- 
nected to IV'o) by any physical operator, and, therefore, 
no spontaneous dimerization is implied in the thermo- 
dynamic limit. Indeed, as shown in Fig. E^b), also in 
geometries with non-zero aspect ratio the p-BCS wave 
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FIG. 4: (a): Difference in the expectation value of the 
nearest-neighbor spin-spin correlations multiplied by the 
number of sites between the two orthogonal states with the 
same quantum numbers predicted by the LSM theorem, (b): 
Finite-size scaling of the dimer order parameter in rectangular 
lattices with finite aspect ratio. The variational parameters 

+ 1. 



are A x z_ y 2 = 2 and A xy = 1 



and L y = L x 



function has a finite dimer correlation length. We have 
therefore shown a clear counter example to the so-called 
Oshikawa conjecture that no spin liquid is possible in 
two or higher dimensions. For the Hamiltonians having 
the wave function as the unique GS, the above 
argumentations are clear cut and conclusive. 

In general, projected wave functions of the type lfT|l. 
with a finite gap Abcs > 0, do not necessarily break 
time reversal, like, for instance, by using a chemical po- 
tential outside the band, d + is symmetry, or 
more simply s symmetry. For all these wave functions 
we expect a finite dimer order parameter in lattices with 
infinite aspect ratio in agreement with the LSM theorem. 

The finite dimer correlation length represents a re- 
markable property of this projected chiral BCS wave 
function. For instance, the conventional short-range 
RVB PjL displays instead power-law dimer correlations 
in 2D, |lj| and therefore gapless features which may de- 
scribe a singular point rather than a 2D spin-liquid phase. 
The fundamental difference between our chiral RVB and 
the short-range one is due to the violation of the Marshall 
sign rule in the former case. This property is observed 
for instance in the ground state of the J\ — J2 model 
in the strongly frustrated region J 2 — 0.5.0 Thus we 
expect that a true spin liquid phase cannot be stabilized 
in particular spin models where the wave function signs 
are not allowed to change as a function of the parame- 
ters of the Hamiltonian. [l9j The complex wave function 
with d-wave symmetry proposed here appears to be a 
reasonable way to open a gap close to a gapless antifer- 
romagnctic phase in a fully translationally invariant spin 
model that allows the violation of the Marshall sign rule, 



i.e. generic frustrated models on a square lattice. The 
chiral order parameter can be indeed small and vanishing 
for A xy — > when the correlation length £ — > 00. A chi- 
ral spin liquid of this form with a reasonable correlation 
length £ ~ 20 provides therefore an almost negligible chi- 
ral order parameter, which is extremely difficult to detect 
numerically and may fit well with the present experimen- 
tal resolutions in High temperature superconductors. |2oT | 
In conclusion we have presented a clear example that 
a spin liquid GS with a gap to all physical excitations, 
though being with single electron per unit cell, can be 
realized without violating the LSM theorem and its gen- 
eralizations in 2D. Q Our results provide a clear support 
to the possibility of a true Mott insulator at zero tem- 
perature, an insulator that cannot be adiabatically con- 
tinued to any band insulator, showing that the effect of 
correlation maybe highly non trivial in 2D systems. 
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